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The large N expansion of matrix-valued field theories was invented more than 
twenty years ago by 't Hooft as a way to treat four- dimensional QCD with gauge 
group SU(N). At N = oo planar diagrams dominate and it was hoped that this fact 
would lead either to analytic results or to a reformulation of QCD as a string theory. 
While the original ideas remain attractive, the program has not yet been successful. 
It was slowly understood that N = oo field theories retain much of the complexity of 
the generic N case. 

Some of this complexity remains even in zero and one-dimensional matrix "field 
theories". In a famous paperB by Brezin, Itzykson, Parisi and Zuber it was demon- 
strated that these so-called matrix models are non-trivial - but still solvable - systems. 
E.g. the zero dimensional model 

VM e~ N Tr [ ^ m2 ~ E ~=i t " Mq \ (1) 

where M is a N x N Hermitian matrix and the t q s are coupling constants parametriz- 
ing a general potential, is solved^ by changing variables to the eigenvalues of the ma- 
trix M and thereby reducing the number of degrees of freedom from N 2 to N. The 
latter proves possible due to the invariance of the above action and measure under 
the group U(N). The model's free energy describes at N = oo a sum over all planar 
diagrams G of spherical topology, where vertices v q of order q (there are #v q of them 
in G) are weighted with a factor t q : 

\o g z ~ z=y: n i i vq ( 2 ) 

G n,eG 

The rather non-trivial combinatorial sum (0) is elegantly calculated by solving ([I]). 
At the timeS the model was considered as a kind of "toy", zero-dimensional QCD 
(retaining the diagrammatic structure but nothing else) to test a new technique. A 
few years later, however, it became clear that the results could be used to obtain the 
solution of a rather interesting physical problem: two-dimensional quantum gravity 
ii'i. Indeed, it was argued that by tuning the couplings t q in an appropriate way 
a continuum limit could be reached at which the planar graphs condense to give a 
continuum path integral over two-dimensional metrics g a b of spherical topology: 

Zoom = / Vg ab e -/ d V^ (3) 

The distance from the critical point in the space of the t q s turns into a continuum cos- 
mological constant /i controlling the area of the surfaces. We also wrote the Einstein 
term, which is however known to be a constant in two dimensions. This approach 
was subsequently worked out and justified by a large number of researchers. 

Thus already the simplest matrix model (|]) contains non-trivial physics. Various 
generalisations were solved and shown to give new physical information; e.g. certain 
simple multi matrix models describe the coupling of c < 1 conformal matter to 2D 
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gravity (||). Such generalisations usually required the development of new techniques 
in order to succeed. Turning this around, each technical advance in large N theory 
usually allows to address a previously inaccessible physical problem. Keeping in mind 
the final goals - like large N QCD or string theories in physical dimensions - a valuable 
strategy is, then, to continue to enlarge our tools and methods. 

An encouraging example for this strategy was given recently in a series of papers&&@. 
We studied a generalisation of ([!]) consisting in the inclusion of an external matrix 
field A in the potential: 

Z = J VM e~ N Tr [ ^ A/2 " s r=a t,(MA)<l (4) 

On a technical level, this model seemed for a long time unsolvable: The external field 
destroys the U(N) invariance of the action. Thus none of the usual methods, whose 
essence consists in reducing the number of degrees of freedom from N 2 to N, appears 
applicable to (§]). We succeeded in nevertheless finding a reduction of the number of 
degrees of freedom by developing the new method of large N character expansions . 
This allows us to address new physical questions. Indeed, it is easy to prove that the 
perturbative expansion in planar graphs of (f|) is given by 

bgz ~ z=y: n tf v *t**«, (5) 

G v*,v q £G 

where 

* = \h Tr A '- (6) 

The sum is, as in eq.(0), taken over all planar graphs G of spherical topology. But 
now we have an extra set of coupling constants t* at our disposal: They assign weights 
to the vertices v* of the dual lattice (there are #v* of them in G) . Thus the t q and t* 
control the coordination numbers of the vertices and the faces of G, respectively. In 
particular, if we set t q — t* = 5 q ^ the only surviving graphs in the ensemble {G} are 
regular square lattices? Thus, the model (§),(0) is capable of interpolating between 
fluctuating random lattices and flat, regular lattices! In the continuum formulation, 
this property would be achieved by adding higher curvature counterterms to the 
action of @: 

Z COQt = I Vg ab e -/^v^ (7) 



Tuning the bare coupling f3$ to zero clearly suppresses any non-flat metric. Our model 
therefore furnishes a precise invariantly regularized definition of 2D higher curvature 
gravity. Before discussing the physics of the latter, let us first explain the steps that 
lead to a solution of the apparently intractable matrix model We will only sketch 



a Of course there is no square lattice of spherical topology; one needs to also add some positive 
curvature defects to be able to close the graph into a sphere; see below. 
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the derivation (and even the results!); a much more detailed discussion can be found 
in the original papers@'!@. 

The idea is to expand the potential in (|]) in terms of the characters of the product 
matrix MA. It is clear that this is possible since the potential is a class function on 
the group (i.e. it only depends on the eigenvalues of the matrix MA). What is less 
obvious (proven in our first work!) is that the expansion coefficients can themselves be 
written as the characters of an auxiliary external matrix B generating the couplings 
t q , just as A generates the couplings t* (see eq.(||)): 

One then has the expansion 

e NTr K~ M MAy = c ^ Xr{b) Xr{ma)) (g) 
R 

with c a numerical constant. The characters are defined by the Weyl formula 

xm(B) = (10) 

where the bi are the eigenvalues of the matrix B, A (b) is the Vandermonde determi- 
nant of the eigenvalues, the set of {h} are a set of ordered, increasing, non- negative 
integers, and the sum over R is the sum over all such sets. The i?'s label represen- 
tations of the group U(N) and the sets of integers {h} are the usual Young tableau 
weights defined by hi = i — 1 + ^boxes in row i (the index i labels the rows in the 
Young tableau, i = 1 corresponding to the lowest row). Note that the restriction on 
the allowed Young tableaux that any row must have at least as many boxes as the row 
below implies that the {hi} are a set of increasing integers: h i+ i > hi. Substituting 
equation (||) into the integral in equation (|j), we can now do the angular integration 
using the key identity 

J (vn) H xr^MQ^A) = d^ 1 X r(M) X r(A), (11) 

where da is the dimension of the representation given, up to a constant, by d^ ~ 
A(h)), and arrive, after performing a Gaussian integral over the eigenvalue degrees 
of freedom, at the expression 

r\ (h e - 

Z = c E n f/,«-ftoV' *W< i4) xw(B), (12) 

{h e ,h°} lLi,jV L i ,b j) 

where c is an immaterial constant. The sum is taken over a subclass of so-called even 
representations. These are defined as possessing an equal number of even weights hi 
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and odd weights h° (since the mentioned Gaussian integration vanishes if the latter 
condition is not satisfied). The formula (0) was originally discovered by Itzykson 
and Di Francesco^ by summing up "fatgraphs" , using purely combinatoric and group 
theoretic arguments. We observe that the matrix model (f|) is thus reformulated as 
a sort of "statistical mechanics model" in Young weight space. The important fact is 
that there are only N weights hf, therefore the reduction to N degrees of freedom is 
achieved. 

The expansion (^) can be further generalised. Consider the matrix model 

Z = J V( p e ~N^ [100+ - £^ =1 g k (^+B)"]^ ( 13 ) 

where is a complex N x N matrix. Introducing a third external matrix field C and 
defining, in analogy with @,(@), 

9k = -— Tr C k , (14) 
yk kN K ' 

one finds by the method outlined above the character expansion 

Z = C E%TY XihM) X{h} (B) X{h }(C), (15) 
{h} ^\ n ) 

where this time the sum extends over all representations. In the special case gu = b~k,2 
the earlier expansion fll2|) is recovered? 

The reduction of the number of degrees of freedom is only a conditio sine qua non. 
One next has to take the large N limit of the expansion. The basic idea is the same 
as for the original model (p]) (see@), with Young weights replacing eigenvalues: The 
weights jfhi are assumed to freeze into a smooth, stationary distribution dh p(h), 
where p(h) is a probability density normalized to one. The details, however, turn out 
to be much more involved. Some simple first examples were worked out in our first 
paper!. An unpleasant feature is that, while the saddlepoint always exists, the support 
of the density p(h) does not necessarily remain on the real axis for completely arbitrary 
couplings t q ,t* complicating the general analysis in a significant way. However, if we 
restrict our attention to models in which the matrices A and B are such that traces 
of all odd powers of A and B are zero the problem does not arise: t% q +i = tjL+i = 0- 
This means that our random surfaces are made from vertices and faces with even 
coordination numbers. Thus it is easiest to consider surfaces made up from squares 
(as opposed to, say, triangles): t* = 5 q ^. A weight t 2q = ^jjTrB 2g is assigned 
whenever 2q squares meet at a vertex (see Fig. 1). 

h For this special case the correspondence between the hermitian model (Q) and the complex model 
(13) was aheady noted previouslytl For a graphical explanation of this correspondence see our 
second workfj. 
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Fig. 1 Flat space and curvature defects. 



This model is clearly capable of describing the transition from flat to random 
graphs. The sum (H) over spherical lattices G4 built from square plaquettes becomes 



G 4 v 2q eG 



(16) 



where vi q are the vertices where 2q plaquettes meet and j^v^q are the numbers of such 
vertices in the given graph G 4 . In order to investigate flT2| ) in the large N limit, one 
attempts to locate the stationary point. This leads to the saddlepoint equation 



2F{h) + i a dti £^L = -hyh 

y ' Jo h-h' 



(17) 



F(h) of the variation of the character 

[ 2 ,}\ . (18) 
where we have also used that the matrix B will satisfy TrB 2q+1 = if we introduce 

N V N 

a y x — 



Here F(h) denotes the large N limit F(hk 
X{h} (B) in eq.(0): 

d 

F(h k )- 1 



dh% 



In 



matrix 62 in terms of which B and the character X{K}{B) are given by 



B 







and X {h}{B) = X{ ^ } (b)x { ^l } (b) sgn[l[(h^ - h°)}. (19) 



i,0 



The variation of the character of the matrix A is easily computed directly because of 
the simple choice t* = 5 q ^. As has been discussed in detail before!, the saddlepoint 
equation (O) actually does not hold on the entire interval [0,a], but only on an 
interval [b, a] with < b < 1 < a: Assuming the equation to hold on [0, a] would 
violate the implicit constraint p(h) < 1 following from the restriction h i+ i > hi. The 
density is in fact exactly saturated at its maximum value p(h) = 1 on the interval 
[0,a\. It is useful to introduce in addition the weight resolvent H(h) as follows: 



H(h) 



h — n' 



(20) 
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We found the weight resolvent H(h) to be very closely related to the standard matrix 
model eigenvalue resolvent (see&). It provides a direct link between the statistical 
distribution of Young weights and the correlators of the model: 

- TiM 2q = - / — h* e^ h \ (21) 
N ql 2m K J 

Here the contour encircles the cut of e H ^ h \ 

The solution of fllTD evidently requires knowing the function F(h). A rather 
general method for its determination has been one of the main technical achievements 
of our work. The method of functionally determining an iV = oo character might 
prove very useful for other applications. We found the following simple result: Further 
introduce the function G(h) as 

G(h) = e H W +F W, (22) 

where, again, the contour encircles the cut of e H( - h \ Its importance stems from the 
fact that it relatesi in a simple way the introduced functions F(h),H(h) and the 
coupling constants ti q '■ 

t 2q = - I ^ G(h)*. (23) 
q J 2m 

It is then easy to deduce, by changing variables from h to G, the expansion 

Q i oo 

*-i = x;^ + £«*G* (24) 

By considering)! the alternative representation of the Weyl character as a determinant 
of Schur polynomials one derives^hat the coefficients a q of the positive powers of G in 
( [24] ) are directly related to correlators of the matrix model dual to (^), i.e. the model 
with t q <-> t*: 

a q = (^Tr (MB) 2 *). (25) 

We have also assumed for the moment that only a finite number Q of couplings are 
non-zero (i.e. t<i q = for q > Q). Furthermore, we were able to show in our second 
worki that ( p4[ ) implies the functional equation 

^W = LiI h -\[G q {h), (26) 

l Q g=l 

where the G q (h) are the first Q branches of the multivalued function G(h) defined 
through ( p4]) which map the point h = oo to G = 0. The resulting picture of the 

c Onc also finds the relation a q = |^g^-lnfx{ii}(^))- As mentioned beforei, rewriting eq.(p4]) 
with this expression for a q suggests a possible relationship to integrable hierarchies of differential 
equations. This conjecture was developed in several conversations with I. Kostov. 
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analytic structure of G(h) is as follows: the couplings ti q determine the number of 
sheets attached to the physical sheet by the cut of e F ^; this number is Q for a 
finite number of non-zero couplings. In turn, the parameters a q determine the sheets 
attached to the physical sheet by the cut of e H< ^ h \ This picture is easily verified (and 
was in fact discovered in this way) for the rather trivial cases where the potential in 
eq.(||) is at most quadratic. The case Q = 2 is shown in Fig. 2. 



cut of e H < h > 




physical sheet 



Fig. 2 Sheet and cut structure of G(h) 



The saddlepoint equation fll7D , together with (|26|) , defines a well-posed Riemann- 
Hilbert problem. It was solveTexacily and in elicit detaii for the case Q = 2, 
where the Riemann-Hilbert problem is succinctly written in the form 

2F(h) + $(h) = -\nh 

21?(h)+H(h) = H- t i), (27) 

The first equation is the saddlepoint equation ([[7]) with I/I{h) denoting the real part of 
H(h) on the cut [b, a] (the righthand cut in Fig. 2). The second equation is ( p6|) with 
f{h) denoting the real part of F{h) on a cut [— oo, c] with c < b (the lefthand cut in 
Fig. 2). This case corresponds to an ensemble of squares being able to meet in groups 
of four (i.e. flat points with weight £ 4 ) or two (i.e. positive curvature points with 
weight £2) (see Fig. 1). We termed the resulting surfaces "almost flat". It turned 
out that all the introduced functions could be found explicitly in terms of elliptic 
functions. E.g. the density satisfying (|17D is given by 

p(h) = -Un ~\<^±,k) with k = J^ 1 (28) 
K \ a — \ a — c 

where sn -1 is an inverse Jacobi elliptic function and K is the complete elliptic integral 
of the first kind (depending on the modulus k). The expressions for the cutpoints 
a, b, c as well as the expressions for H(h), F(h), G(h) and the physical correlators a q 
(see (^)) were found as well.] 
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The resulting surfaces are very different from pure gravity. The spherical partition 
function consists of flat cylinders pinched shut at the two ends, resulting in 
precisely four defects. Very short, highly twisted cylinders dominate the continuum 
limit. By exactly calculating the correlators a q (see eq.(|25|)) we also analytically solved 
the combinatorial problem of surfaces with a single negative curvature insertion of 
arbitrary degree balanced by a gas of positive defects. 

To analyse the problem of the transition from flat to random lattices, we merely 
need to perturb our almost flat lattices by any operator containing negative curvature. 
This physical observation allows us to extend the Q = 2 solution in a simple way. 
Choosing the weights to be 



to — v^A t. td — A, tfi 



A§ ^ 
A t ' 



.. t 



2q 



V t 



the expansion (24) becomes 
h 



1 = ^ + 

G G (G — e) 



+ positive powers of G, 



(29) 



(30) 



where e = |j/3 2 . The derivation of the functional equation ( p6[ ) is easily modified 
(note that we now have an infinite number Q = oo of weights) to give 



H(h) 



h 



tto — La 



Gi(h) G 2 (h). 



(31) 



The essential point is that one keeps the property that only one other sheet 
is attached to the physical sheet G\(h) by the cut of e F ^ (see again Fig. 2). The 
only difference is that the semi-infinite cut [— oo,c] becomes a finite cut [d,c]. This 
results in modifying the Riemann-Hilbert problem ( p7|) to 



(32) 



2F(h)+fi(h)= -\nh 
2f(h)+H(h) = hi(^^) 

It may still be explicitly solved in terms of elliptic functions; e.g. the density (|28| ) is 
generalised to 



p(h) 



u 
K 



In 



7T 



Oij^ju -iv),q) 
[u + iv),q) 



0i( 



2 A' 



with 



and k 



(a — b)(c — d) 



(33) 



(34) 



\ (a-c)(b-dy 

where K and K' are the complete elliptic integrals of the first kind with respective 
moduli k and k' = a/1 — k 2 and u and v are given by 



u 



sn 



(a-h)(b-d) 
\ (a-b)(h-dy 



k) and v = sn 1 ( 



I a — c 
a — d' 



k'). 



(35) 
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Again, the cutpoints a, b, c, d can be explicitly obtained as functions of the couplings 
X,(3,t. 

Let us also mention that the weights can be further generalised while keeping 
the quadratic (two-sheeted) structure of the function G{h). The idea is to shift away 
the simple pole at G = in the expansion (|3"0D, i.e. 

h — 1 = — h — h positive powers of G. (36) 

G — €± G — €2 

Here ci and Co are two constants given by ci = t262 ~ ti and c 2 = u ~ t2ei . It leads to 
the Riemann-Hilbert problem 

2F(h)+l/l(h)= -ln/i 

2f (/i) + = ln( (£l+e f 2 - t4 + ei e 2 V), (37) 

The weights generalizing (|29|) by one extra parameter can be found explicitly by 
expanding (^TJ) in inverse powers of G. We have not worked out this further explicitly 
solvable case in detail. Analysing it would furnish an interesting universality check 
of our result. 



Having at hand the explicit solution of the model for the weights (29), we are in 
a position to analyse the problem of discrete 2D R 2 gravity. With these weights it is 
easy to prove, using Euler's theorem, that the partition sum ( |16D becomes 



Z(t,\,P) =t A £ X A (3 2{#V2 ~ 4 \ (38) 

G 4 

where A is the number of plaquettes of the graph G and j^v 2 the number of positive 
curvature defects. Note that the latter are balanced by a gas of negative curvature 
defects, whose individual probabilities are given in (|29|). 

We expect this model to describe pure gravity in a sufficiently large interval of /3, 
after tuning the bare cosmological constant A (controlling the number of plaquettes) 
to some critical value A c (/3). On the other hand, for A fixed and (3 = we entirely 
suppress curvature defects except for the four positive defects needed to close the 
regular lattice into a sphere. It is thus clear that (3 is the precise lattice analog of 
the bare curvature coupling (3 in the continuum path integral ([7p. The phase (3 = 
of "almost flat" lattices - very different from pure gravity - studied in detail in the 
second papeiH was discussed above. 

Let us now summarize the main physical results following from the exact solution 
(for general A and (3) of this model: 

1. A long debated question was whether models of the present type undergo a 
"flattening" phase transition at a finite, non-zero critical value of f3 = f3 c . The weak 
coupling region (3 > (3 C would then correspond to the standard phase of pure gravity 
while a putative novel "smooth" phase of gravity might exist either at (3 = f3 c or in 
the entire interval < (3 < f3 c . This would constitute an existence proof of continuum 
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2D R 2 gravity. We found analysing the exact solution, to the contrary, that there is 
no "flattening" phase transition at non-zero (3. For any given (3 we find the powerlike 
scaling of standard pure gravity on large scales. This means that no matter how flat 
the system is on small scales (of the order of it destabilizes in the infrared into 

the familiar ensemble of highly fractal "baby-universes" . 

2. The dependence of the partition sum (|38D on f3 and the lattice cosmological 
constant A in the vicinity of the flat phase /3 ~ and close to A ~ A c is given by 
a simple, (presumably) universal scaling function f(x) (defined through Z(t,\,/3) = 
f(x)) reflecting the transition from flat space to pure gravity: 

/(a;) = x 6 --x 4 + — x 2 - — -x(x 2 -l)5, (39) 
2 8 16 v ; v ; 

where the scaling variable x is given, to leading order, by 



ML 

15/3 2 



*=^i^. (40) 



7T (3 



We can distinguish the following features: 



(a) There is a degree | singularity at x = 1, correctly^ reproducing the universal 



string susceptibility exponent j s = — \ of pure gravityoH. In view of eq. (flCf) , the 
critical value of the lattice cosmological constant A is therefore given to leading order 
by A c = 1 — -^(3 + 0{(3 2 ). Therefore (see (|38|)), the characteristic growth of the 
random surfaces as a function of the lattice area A (= number of plaquettes) is given 
by 

Z(t,A,(3) ~ t^ e ^ pA A-i. (41) 

For any non-zero (3 we do have exponential growth of the number of surfaces, but one 
has to go to larger and larger scales (i.e. use more and more plaquettes) to be able 
to take the continuum limit. If (3 is exactly zero there is no longer any exponential 

g 

growth and no pure gravity continuum limit is possible. The prefactor (3~z is found 
in the exact calculation in section 5; we are not sure whether it is universal. 

(b) We further see that taking /3 — > before the limit A — > A c corresponds to the 
limit x — > oo. In this limit one finds f(x) ~ jJs ^ + C(^) ; that is, the characteristic 
critical behavior of 2D gravity "silently" disappears and we recover a power series in 
- corresponding, in view of (0), to a perturbative expansion in lattice defects (3. In 
this limit the characteristic growth of surfaces as a function of area A is 

Z{t,A,P) ~ t 4 ( A + 0{(3 2 A 3 ) ). (42) 

The leading order corresponds precisely to the almost flat lattices (with exactly four 
positive defects) studied in our second work!. The corrections are interpreted as 
insertions of negative defects, balanced by further positive defects. The typical shape 
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of the surfaces in this limit is a generalisation of the one we found for "almost flat" 
graphs: Long, filamentary cylinders growing out from every negative curvature defect. 

(c) It is easy to prove that the scaling function fl39|) is the simplest possible function 
with the limiting properties discussed in (a) and (b). 

The above results might be interpreted in terms of a continuum model of quantized 
curvature defects, in which the localised defects move around like particles in a gas 
on a flat background space. The deficit angle, A9, of a defect can take on the values 
AO = n,0 and — tt. A positive curvature defect is surrounded by a conical geometry, 
whereas a negative curvature defect corresponds to a saddle-type insertion (see Fig. 1). 
The higher order negative curvature defects (— 2tt and higher) would not be expected 
to play a role in this limit (the entropy from moving two low order defects around 
would completely dominate that from a single higher order defect). The coupling (3 
can be interpreted as a fugacity controlling the number of defects. The flat space limit 
(3 —>■ consists of four defects of degree 7r moving around with respect to one another. 
Varying the fugacity, (3, allows one to smoothly interpolate between flat space, (|42|) 
(with four defects), and pure gravity (|4l| ) (with an infinite number of defects). 

One might also attempt to develop this picture directly in the continuum. One 
could start with the conformal metric of a flat surface with localised curvature defects. 
It can be represented locally as g a b = 5 a b e^ z > with 

M 

cp(z) = £ R j ln(z - z 3 )\ (43) 

where Rj = —1,1. Symbolically, the partition function might be written as 

Z(ji,) = $ M I d[z!,...,z M ] e~^ d2z V det s«. (44) 

Here we introduced the fugacity of curvature defects (3 instead of the explicit i? 2 -term 
in the action. It serves the same purpose: for (3 —>■ we arrive at the completely flat 
metric, whereas for j3 ~ 1 the system should show the behaviour of pure quantum 
gravity, at least in the infrared domain. We retained the notation f3 to denote the 
parameter playing a role similar to the R 2 coupling in the above discrete model. 

This formulation resembles a little bit the two-dimensional Coulomb gas problem. 
However, the measure of integration d[z\, Zm) of the positions of the curvature 
defects is a complicated object: it should take into account the topology of the surface 
and the existence of zero modes (the action does not depend on some directions in 
the space of the Zj). It would be very interesting to make this direct continuum 
formulation more precise. 

Another interesting issue is the role of exponential corrections appearing due to 
the structure of elliptic functions. In fact all physical quantities, such as Z(t,X,/3) 
(see eq . (|3~8D ) , contain exponentially small terms in the limit A — > A c and (3 — > 0, thus 
leading to an essential singularity at /3 = 0, A = 1. One can obtain the first correction 
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of this type in e.g. the free energy f(/3) per unit area in the thermodynamical limit 
A A,.: 



z(t,A,(3) ~ t mA A ~\ with / ( ^ ) = _L^r (1 + ... ) + e -^ (4 + ...)] (45) 

where /(/?) = lim^oo i \nZ(t, A, (3) = A c and the dots denote terms of order /3 3 and 
higher. 

These exponential terms are likely to be lattice artifacts. They emerge even in the 
simplest calculation for the flat closed quadrangulation with four positive curvature 
defects, where they appear as discrete corrections to the approximation of elliptic 
sums by integrals close to the continuum limit. 

On the other hand, formula (^) corresponds to the critical free energy as a func- 
tion of the curvature fugacity (3 (i.e. we have already taken the continuum limit). It 
is possible that the exponential terms might be corrections relevant for the statistical 
mechanics of random lattices at long distances (of order -s) rather than for continuous 
2D-gravity. 

In conclusion, we have tried to emphasize two points in this presentation: 

1. Technical advances in large N group theory are possible. Apparently hopelessly 
difficult problems like the model @ become tractable when one changes the technique. 

Our approach raises a number of interesting questions. First of all, the method 
we have employed seems rather indirect and is definitely very involved. Is there a 
hidden simplicity we have missed so far? Are there alternative ways to obtain our 
results? Secondly, one should ask whether the method could be adapted to other 
theories describing new interesting physics; in particular gauge theories and strings 
in physical dimensions. 

2. Each such advance leads to the capability to address new physics questions. The 
present advance enables us to treat the problem of two-dimensional higher curvature 
gravity. 

We have presented for the first time an exactly solvable model interpolating be- 
tween flat space and 2D quantum gravity. This addresses a long-standing open prob- 
lem. We would like to stress, however, that we have by no means demonstrated yet 
the universality of our result. Could it be that finetuning the weights t q leads to new 
phases of gravity? Here we certainly do not mean to rederive the usual multicritical 
phases of the one- matrix model (|I|). We rather envision tuning the weights in a subtle 
way so as to reach a phase of smooth 2D gravity. Certainly the analytic complexity of 
the result could be a hint that much more waits to be discovered. (Compare e.g. the 



densities (|28|) , (|33|) of the model ([4|) under study with the rather trivial algebraic den- 
sities obtained in the absence of the external field.) Before addressing this issue, we 
repeat, one first has to simplify the method enough to allow for a deeper insight into 
the analytic structure of the solution for general couplings t q . 
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